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Abstract 

In the framework of QCD at finite temperature we have obtained dissipative terms for the 
effective potential between q and q which would partly explain the J/il) suppression in the Quark 
Gluon Plasma (QGP). The derivation of the dissipative potential for QGP is presented and the 
case for Hadron Matter (HM) is briefly discussed. The suppression effects are estimated based 
on simple approximations. 
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1. Introduction 


It was supposed that a suppression of J/dl production in Relativistic Heavy Ion Collisions 
(RHIC) could be a clear signal for formation of a new matter phase, e.g. the Quark Gluon 
Plasma (QGP) |l[], later the phenomenon was observed by the NA38 collaboration [Q. However, 
it is too early to celebrate the discovery of QGP, because rigorous studies indicate that in the 
Hadron Matter (HM) phase, the suppression phenomenon also exists. Thus many authors study 
the mechanisms which cause the suppression of J/T production in the QGP and HM phases 
separately |l| 0. They investigate how the mass of J/T shifts with temperature and density of 
the QGP and HM medium, especially the dependence of the quark and gluon condensates on 
temperature and density are introduced to explain the mass shifts and then the suppression of 
J/T production. 

The RHIC experiments provide an atmosphere of hot and dense matter state, no matter it is 
in QGP or HM phases. The two situations are quite different in many aspects and an important 
research subject is to look for one or several clear signals to conhrm or negate formation of QGP. 
The J/T suppression is one of the candidates, but obviously, observation of such suppression is 
not enough for drawing a dehnite conclusion, one should investigate the details of the mechanisms 
which result in the suppression in QGP and HM. If one starts with a fundamental principle, 
i.e. QCD, instead of phenomenologically obtained potentials, he may expect to derive different 
potential forms for QGP and HM, and it just is the task of this work. 

The effective potential between quark-antiquark or quark-quark can be written as 

V (r) =+ ar. (1) 

r 

In the regular zero-temperature theory Oe// = for meson and 2as/3 for baryon where 

as = Qs/'i'n' and Qs is the QCD coupling constant. At hnite temperature, the conhnement 
constant becomes temperature-dependent [^], 

^ dec 

where T^ec is the deconhnement temperature and 6 is an uncertain parameter. Above T^ec, the 
linear conhnement potential disappears, but it does not mean that the bound state dissolves. 
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because a binding energy can also be provided by the Coulomb-type potential. By analysis of 
the medium state, a Debye screening mechanism is suggested as 

V (r) = (3) 
r rn 

where is the Debye screening length. It is a function of temperature and may be obtained 
by models. The dissolution condition can be written as 

rn < OMtb = 1.68/{aeffrUc). 

Seeking for other possible mechanisms which may cause the J/T suppression, we turn to 
re-study the situation. We will start with the QCD theory and its low energy phenomena at 
hnite temperature and density. 

Let us briefly retrospect how one derives the Coulomb-type potential from the held theory. 
Considering a t-channel scattering between two quarks or quark-antiquark, the amplitude in 
momentum space is 

M = ui{pi)j^^U2{p2)G>"''{k^)u3{p3)'yuUi{pi), ( 4 ) 

where is the full gluon propagator and for g — g interaction and M 4 become and v^. 

There is also an extra s-channel annihilation diagram for q — q, where U2 and M3 turn to V2 and 
M 3 . Setting fco = 0, a Fourier transformation of the amplitude gives rise to the effective potential 
in conhguration space. Generally speaking, in t—channel the momentum transfer is space-like, 
i.e. k^ < 0 . 

It is known that the photon and gluon propagators cannot provide an imaginary part, even 
at hnite temperature and density in the common sense. The reason is that an absorptive part 
of the photon and gluon propagators corresponds to on-shell photon or gluon (gauge bosons), 
i.e. = 0, it indeed manifests an emission or absorption of a real photon or gluon from the 
quark (antiquark). At s-channel, the intermediate gauge bosons have momentum k^ = s = 
iPq + Vq)'^ > > 0, thus the = 0 condition cannot be satished, unless via loops. Whereas 

at t-channel the momentum transfer carried by the gauge boson is = {pq — p^)^ < 0, so it 
seems that an imaginary part can appear. But it is not true, because only at the boundary point 
iPq ~ Pq)^ = 0 , it means that if at both sides of the interaction, the quark is on its mass shell 
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Eq = Eg and Pq = i.e. the 6 = 0 forward scattering. For the Fonrier transformation the 
integration should be restrained in a momentum conservation allowed region, so for the exact 
on-shell situation, the integration region is zero (see below in the context for details), namely 
lo l^l/(^) would result in a null contribution to the absorptive part. By contraries, for the real 
part of the propagator which gives rise to the regular potential, = 0 condition is not imposed, 
so the integration region does not have any restriction (see below). That is why in the regular 
theories, photon and gluon propagators do not contribute an imaginary (absorptive) part, but 
only a real (dispersive) part. 

However, in a bound state under consideration, the situation is different. If the heavy quark 
or antiquark (c, c) is on-shell, it can be described as 


p^ = rricV^, (5) 

where is the four-velocity, but in a bound state there is an integration region [Q 

p^ = rricV^ + k^, ( 6 ) 

where is the ’’residue” momentum, and \k\ ~ fco < ^qcd- Thus the integration region allowed 
by the energy-momentum conservation is no longer zero, but can be from zero to Aqcd (or a 
smaller value than Aqcd)- If can indeed contribute an absorptive part which finally results in a 
dissipative term in the effective potential. In the next section we will present the details. 

Assuming the temperature and density are above Tdec and pdeo the QGP phase is reached. 
Without the linear conhnement term whose source is purely non-perturbative and obscure so 
far, we only need to take into account the one-gluon-exchange diagram which results in the 
interaction between quarks and quark-antiquark. It is well known that at T = 0, the leading 
term of the one-gluon-exchange contribution is the familiar Coulomb-type potential —aeff/r, 
but as the non-perturbative QCD effects which are characterized by the non-vanishing quark 
and gluon condensates are taken into account, the potential form is modified and corrections 
related to l/mg {rric or mb) emerge 00]. When the temperature is non-zero and the density is 
above the regular one, more extra contributions appear. In this work, ignoring all spin and spin- 
orbit dependent terms and under some simple approximations, we obtain a new Coulomb-type 
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potential as 


V{r) = + i(a + asb)], 


(7) 


where a and b are functions of temperature and density. This new term turns the J/T charmo- 
nium into a dissipative system and a dissolution is expected. 

In the next section, we derive eq.(|^ in the framework of QCD at hnite temperature and in 
Sec.Ill, we discuss its signihcance and the situation for HM phase. 


II. Formulation. 

In the QGP atmosphere the propagator of fermion at hnite temperature and density in 
momentum space can be written as B 


where 


Spik) = (^ + m)[ 


fpix) = 


^2 _ ^2 


9{x) 


27i6{k‘^ — m?)fp{k ■ u)], 


9{—x) 


( 8 ) 


(9) 


-|- 1 + 1 ' 

In the expression (3 = 1/T, u is the four-velocity of the medium, generally u = (1, 0) in the 
laboratory frame, so k ■ u = ko, fi is the chemical potential and is related to the density of the 
medium. Whereas the propagator of gluon reads 


where 


D^u{k) = [-p + 27i6{k‘^)fB{k ■ u)]{-g^^ + 


fsik ■ u) = 


g/3^0 — 

It is noted that we are working in the real-time scheme 


( 10 ) 


( 11 ) 


As aforementioned, if the charm-quark in J/T bound state is strictly on mass shell, the t- 
channel scattering demands < 0, so that the fermion propagator does not have the temperature- 
dependent term because the intermediate fermion of non-zero mass cannot be an on-shell real 
particle. However, the soft gluon interaction can make the charm quark deviate from its mass 
shell, namely the k‘^ can be equal to or greater than zero in eq.(^ and it represents an emission 
or absorption of real photon or gluon, (/c^ = 0 and \k\ ~ Aqcd), moreover a small time-like 
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momentum transfer > 0 is also reasonable, then the intermediate fermion is an on-shell 

real particle. 

Omitting the loop corrections, the lowest order Feynman diagrams are depicted in Fig.l. In 
general, the contribution can be decomposed into a form 

/(/x,T) = /i(r) +^V2(T) + ^V3(T,/i), (12) 


where /i(T),/ 2 (T) and are from one gluon exchange, gluon condensate and quark con¬ 

densate diagrams respectively shown in Fig. 1. It is noted that only /s depends on the density 
of QGP, but not /i and f 2 - Both g‘^f 2 (T) and belong to higher order expansion, but 

while discussing the density dependence, fsiT, g) becomes the leading term, so we drop out / 2 (F) 
but keep f^{T,g) in our calculations. The second term with the 5—functions in the propagators 
correspond to the on-mass-shell intermediate particle, therefore turns to be a real particle, so it 
can feel the influence from the surrounding atmosphere, namely the density and temperature. 


The traditional way for obtaining the potential is only to deal with the off-shell part of the 
propagator, namely to set fco = 0 and carry out a three-dimensional Fourier transformation of 
the two quark scattering amplitude. All the details including the loop corrections and non- 
perturbative effects are presented in literatures ii 00 , instead this work is only focused 
on the temperature and density-dependent parts which would contribute additional modihcation 
terms to the potential, we denote them as V^ir) and {r) corresponding to Fig.l (a) and (b) 
respectively. 

From Fig.l (a), one has 

V^(r) = 

G\ ) ^ \ 


-IQ'Ki 


j/® 


as){-2v:) 


d^k 


6{k^)e 


2\ -ik-(x-y) _ 




( 13 ) 


where the factor coming as a common factor due to the color singlet condition of 


hadrons and it makes the leading Coulomb term be 


4q:s 

3r 


as required. It is also noted that here 


one does not need to set fco = 0 as in ref.[IT| and the hnal integration region for \k\ is from 0 
to Aqcd as discussed in section I. The factor jA — guarantees the dimension of Voir) right and 


j/» 


comparing with the traditional way for deriving the potential, there is an extra integral over ko, 


so a factor 4 


M 


j/9 


is needed. This is equivalent to multiply a time factor r (r ■ ~ 1). 
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Taking the spontaneous requirement xq = yo (potential means an instantaneous interaction) 
a straightforward calculation gives 


V^ir) = 


—Idvrz 




E 


3 "" (27r)2 rMj/^ 

[r(l — e“”^^cosAr) — sin Ar], 


X 


(14) 


where A < Aqcd is a parameter and in the expression we deliberately pull out the 1/r factor 
(see below). It is a series which absolutely converges for any hnite r and /?. If r is not very 
large (a few tenths of fm in our case), one only needs to take hrst a few terms for numerical 
computations. One can notice that if A ~ 0, i.e. for quarks are exactly on mass-shell, V^(r) 
vanishes and at T —>■ 0, it is also zero, this is consistent with our common knowledge. To see its 
meaning and avoiding tedious calculation, setting A —> cxo we approximate eq. (JT^) to an integral 
as 


E 


- 1 - Ji 


r 1 TT j3 

dx „ -r = —(- arctan—). 


Thus the Vq (r) is recast as 


t^T/ n —IQm 1 1 TT P 

Vr v) = - Cis-, —rxT;-(- arctan—). 

^ 3 (27r)2/?r^2 


Below we will evaluate its contribution in terms of the series solution (^). 
Similarly, for the quark contribution shown in Fig.l (b), one has 


V,^{r) = 


-327ri .^<ipqipq>T 1.1 

r 
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cos yr 


r 




n=l 


^ 2/^2 _j_ ^2 

[2r cos yr — [r cos Ar -\- n(3 sin Ar) (1 ) ], 


X 


( 15 ) 


for A > /i and 


Vf(r) 


—327ri ‘Ipqlpq >T 1 .1 — COS/ir 


-a 


-E(-i)’ 


j—nfiii 


9 * niq r'' r ^ n‘^f3'^ + r‘^ 

[r(2 - (e"^^ - e-"^^) cos Ar) -1- n(3smAr{e^^^ - 


X 


( 16 ) 


for A < p, it is obvious that as A —0, Vq (r) vanishes. While deriving eq. (p!5[) , we approximate 


E = + m2 Ri |/c| in the integral due to the smallness of niq. In the expression < 'tpq'ipq >t is 

the quark condensate at hnite temperature T and its signihcance will be discussed in the next 
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section. 


III. Discussions. 


The additional potential terms induced by QCD at finite temperature and density are imag¬ 
inary and tend to zero as the temperature and density approach to zero. It makes sense because 
if there is no hot and dense atmosphere for the gluon and quarks, the additional terms do not 
exist at all, so this scenario would not affect the regular lifetime (about 7.5 ± 0.4 x 10“^^ sec.) 
of the J/T which mainly is determined by the s-channel annihilation process of qq into three 
gluons. When T is very high (a few hundred MeV) and p S> po, the damping factor becomes 
substantial and it is the QGP situation, (see below) 


The temperature and density effects turn the Hamiltonian of J/T into complex, and the new 
Hamiltonian describes a dissipative quantum system. Because it is dissipative, the system would 
dissolve after a certain time, for example, into DD or DgDg etc. Dissolution to DD or DgDg 
needs to absorb energies from atmosphere because < 2Md, in QGP case, it is possible, 

whereas in the vacuum circumstance, due to the kinematic constraint of the final state phase 
space, even there were a dissipative term, the dissolution would not occur. 


From eqs.([I^, [T^) one can notice that the coefficients of i in VQ{r) and (r) are always 
negative. This is a very important point because it makes the system dissipative. 


Since V)J(r) and {r) are very complicated functions of r, it would be extremely difficult to 
solve the Schrodinger equation with such a potential. Therefore to see the physical significance, 
we would choose a simple but reasonable approximation. 


(a) Estimation of the dissipation effects. 

(i) Looking at the additional Hamiltonian, one can note that they can be written as f/q(G)(^)“ 

where fq(G){f’) are complicated functions of r and given in eqs.(^, ^). Our approximation is 
to treat fq(G){.'^') as average values instead of functions of r, namely, we approximate fq(G){.‘>^) as 
fq{G)if) where f is the average radius of J/T. Matsui suggests 0.2 < < 0.5 fm. 

(ii) As (t{T) disappears above the deconfinement temperature T^eo the stationary Schrodinger 
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equation becomes 


Aa. 


[1 + i{a + as6)]0 = E(() 


2mred 3 r 

where 0 is the wavefunction of J/A/, nired is the reduced mass which equals mc/2, and 

1 


1 OO 

= rE 


VT ?7.2/?2 + 


X [r(l — e cos Ar) — n/ 9 e "^^siuAr] 


A—>oo 1 1 , TT /3. 

—> — — I - arctan—) 

71(3^2 f’ 


b = 


b' = 


„ 8 < 'i’qjjq > A COS/ir ^ Hxn j- 

^ 3 f ?7,2/?2 + f2 


X 


q=u,d,s 


[2r cos/ir — (r cos Ar + U/d sin Ar)e "^^(1 + e^’')], for A >/r; 
8 <^|Jq^pg> 1-cos fif ^ 

^ 3 f 


X 


q=u,d,s 


( 17 ) 


(18) 


(19) 


[f(2-(e"'’'' + e-”'’'')cosAf)+n;3siiiAf(e"'’''-e-"'’'')l for A </i. (20) 


The eigenenergy of eq. m is simple and reads 

(EjUf = ( 21 ) 

where {as)eff = «s[l + i{cL — Usb)]. For the ground state j = 1. 

(hi) The evolution of the quantum system can be expressed as 

= [1 - (a + asb)‘^]t) ■ + asb)alt). (22) 

y y 

There exists a damping factor. In the calculations, we ignore the temperature dependence of 

H- 

(b) Numerical evaluations. 

(i) If the QCD expansion (including the condensates) converges, \bas\ should be smaller than 
a. Typically, ~ 0.3 in the potential model, the term agb can compete with the term a. b is 
proportional to the quark condensate at hnite temperature < 4’q4’q >t which decreases with the 
increase of temperature. However, we will argue in the following that the ratio where 

rriq 

mq is the constituent quark mass is almost independent of temperature. 


















The constituet quark mass is defined as the pole of the quark propagator 


E(p^ = ml) = m„ 


(23) 


where is the dynamical mass related to quark condensate |l6 






< > ««(/)• 


(24) 


Hence the ratio depends only on agim^a) which dependence on temperature is ignored. 

TTiq H 

Therefore, the parameter b is almost temperature independent. 


(ii) Thus we have 

|0j/^(^)r = = 0)pe^ "0 \ 

A typical time factor is about tq ~ 7 x 10“^^ sec., for 0.2 < f < 0.5 fm and here we use f = 0.4 
fm. Even though this numerical value cannot be taken very seriously, the order of magnitude is 
reasonable. 


The size of the collision region is about 10 ~ 100 fm, a particle produced at the center of the 
region needs 1 x 10“^^ ~ 10“^^ sec. to travel to the boundary, so Tq is of the same order as the 
traverse time and definitely results in a suppression of J/T production (in fact a bulk decay) 
observed in experiments. 


In fact, in such a small time interval, according to the regular theory c and c hardly combine 
into a bound state, so it is equivalent to the effective screening. 

(iii) It is worth noting that only b depends on the density via the chemical potential /i. From 


eq. ([l5|) , numerical results show that as |p| gets larger the contribution of {r) becomes more 
competitive to Vq^t), but by the common sense, |6 q;s| < a. 

(c) For the hadron matter phase, gluons and quarks do not directly feel the temperature and 
density of the hadron medium, therefore the gluon propagator cannot be influenced by the 
medium. In this case the temperature and density effects can appear in two ways. One is that 
the quark and gluon condensates are modified in the medium but different from that in QGP, 
pT[| , while another way is via loops in s-channel. 

Our results show that the temperature and density effects can cause a suppression of J/T 
production, but the mechanism is different from the traditional Debye screening effect, i.e. it 
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contributes a dissipative part to the potential. This additional term makes the quantum system 
dissolve by a time scale 7 x 10“^^ sec. This mechanism can only exist in QGP phase but not in 
the HM phase. 
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Figure Caption 


Fig.1.The Feynman diagrams of 


t-channel c — c scattering where the charm-quark may deviate 


from its mass shell a bit. 
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